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ABSTRACT 


This  report  includes  a  theoretical  discussion  of  energy  loss  rates 
and  transport  equations  in  thermionic  converters,  as  well  as  spectrosi opic 
measurements  on  a  cesium  diode  and  a  determination  of  the  electron  impact 
ionization  cross-section  in  cesium. 

Theoretically,  effective  coefficients  for  ionization  and  radiative 
energy  loss  have  been  computed  for  different  values  of  electron  temperature 
and  density.  These  coefficients  and  the  neutral  density  yield  the  rates  for 
net  ion  production  and  radiative  energy  loss  from  all  sources  except  the 
resonance  lines.  Energy  loss  from  resonance  lines  and  the  ion  production 
cost  are  considered  separately.  The  following  general  conclusion:,  can  be 
drawn  from  the  consistent  benavior  of  the  results:  (1)  the  area  within  0.4 
mm  of  the  emitter  is  a  source  of  ion-electron  pairs,  the  rest  of  the  plasma 
being  a  sink;  (2)  inelastic  energy  losses  decrease  monotonically  from  the 
emitter  to  the  collector;  and  (3)  these  effects  are  relatively  independent  of 
cesium  pressure,  but  are  generally  proportional  to  output  current  density. 

furthermore,  the  differed  ial  equations  governing  the  transport  of 
particles,  momentum,  and  energy  in  i  o  plasma  between  the  electrodes  of  a 
th< — ionic  converter  have  been  formulated  in  a  more  complete  fashion  than  has 
been  done  previously.  Specific  advances  are  the  inclusion  of  more  exact  coefficients 
for  electron  transport  and  for  ion  production  and  radiative  energy  loss  in 
the  plasma  volume. 

In  the  diode  experiments,  the  electron  temperature  Tg  was  deter¬ 
mined  from  measurement  of  the  continuum  emitted  by  radiative  recombination 
to  the  6P  states  of  cesium?  and  the  electron  density  Ng  was  determined  from 
Stark  broadening  of  the  fundamental  series  lines.  By  transversing  the  mono¬ 
chromator  with  respect  to  the  discharge,  spatial  measurements  of  Ng  and  Tg 
were  made.  Of  particular  interest  is  observation  of  forbidden  linel.  The  6 
7269  A  line  (6G  -  has  a  very  definite  structure  in  the  middle  and  sharp 

fall-offs  in  the  wings;  its  shape  depends  strongly  on  electron  density. 

The  total  ionization  due  'co  the  passage  of  an  electron  beam  in  cesium 
vapor  has  been  measured  with  a  Tate  and  Smith-type  apparatus.  The  retarding 
potential  difference  method  was  used  in  the  electron  gun  to  obtain  energy 
resolution  better  than  0,1  eV,  The  density  of  Cs  atoms  was  determined  from 
the  Taylor  and  Langmuir  tables  and  also  measured  with  a  surface  ionization 
detector.  The  two  methods  give  the  same  density  to  within  +  3%  when  the 
apparatus  is  t  thermodynamic  equilibrium  with  the  Cs  reservoir.  The  cross 
section  for  production  of  Cs+  ions  has  been  determined  from  threshold  to  100  eV. 

For  energies  above  50  eV  the  results  agree  closely0with  those  of  McFarland 
and  Kinney^  At  threshold  we  found  a  slope  of  2.7  A  ^/eV?  as  compared  to  the 
value  2.2  r.  /eV  determined  by  Heil  and  Scott. 
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1.0  INTRODUCTION 

This  is  the  Second  Annual  and  also  the  Final  Report  on  Contract 
Nonr  5154(00),  covering  the  period  from  1  November  1966  to  31  October  1967. 

The  objective  of  the  work  performed  under  this  contract  has  been  to  determine 
the  properties  of  the  cesium  plasma  in  the  ignited  (or  arc)  mode  of  operation 
of  a  thermionic  converter. 

The  experimental  part  of  the  program  has  involved  spectroscopic 
measurements  to  determine  electron  temperatures  and  densities  in  the  inter¬ 
electrode  plasma  of  a  converter  and  to  explore  features  of  the  spectrum  for 
use  in  refined  diagnostic  techniques.  Ionization  cross  sections  have  also 
been  measured  in  a  separate  experiment,  with  particular  attention  on  the 
threshold  region.  Theoretical  work  has  developed  transport  equations  that 
properly  account  for  interaction  between  the  three  components  (ions,  electrons 
and  neutral  atoms)  of  the  cesium  plasma.  These  equations  allow  the  prediction 
of  variations  of  electron  density  and  temperature  when  initial  conditions  are 
given. 

The  present  work  has  contributed  to  the  growing  knowledge  of  plasma 
behavior  in  converters.  The  use  of  spectroscopy  in  thermionics  research, 
initiated  at  Los  Alamos,  has  been  a  major  factor  in  achieving  our  present 
understanding.  There  is  still  much  to  be  done  to  refine  spectroscopic  tech¬ 
niques.  Forbidden-line  intensities  offer  further  information  about  plasma 
conditions  and  have  not  yet  been  utilized.  The  data  on  forbidden  lines  pre¬ 
sented  in  this  report  are  the  first  quantitative  results  that  we  know  of  for 
cesium. 

Transport  of  electrons,  energy,  and  momentum  across  the  plasma  is 
fundamental  to  converter  operation.  Inelastic  and  elastic  rates  for  energy 
exchange  between  electrons,  neutrals  and  ions  are  now  quantitative  enough 
to  attempt  solution  of  the  relevant  equations.  This  program  has  contributed 
to  our  knowledge  of  these  rates,  and  in  this  report  the  transport  equations 
have  been  formulated  for  computer  solution. 

Satisfactory  understanding  of  plasma  effects  in  thermionics  is  the 
practical  goal  of  much  continuing  research.  Towards  that  goal  this  work 
has  been  dedicated. 
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2.0 


THEORY 

The  theoretical  program  has  been  directed  toward  describing  the 
transport  of  particles  and  energy  across  the  plasma.  This  is  a  part  of  the 
more  general  problem  of  predicting  current-voltage  characteristics  of  a  diode. 

The  first  year's  program  defined  a  model  for  the  interaction  of 
the  ions,  electrons,  and  atoms  in  the  plasma  and  determined  effective  ion¬ 
ization  and  recombination  rates.  The  present  (second)  year's  program  has  gone 
considerably  further.  Rates  for  energy  loss  of  the  electrons  have  been  cal¬ 
culated  and  the  equations  for  transport  of  particles,  momentum,  and  energy 
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have  been  formulated  for  typical  "high  density"  (N  =“  10  cm  )  diode  con¬ 
ditions.  A  rather  complete  description  of  the  plasma  in  the  high  density 
diode  is  now  available. 

2.1  Electron  Kinetic  Energy  Loss  Rates 

Consider  a  steady  state  plasma  at  a  given  electron  density,  electron 
temperature,  and  total  pressure.  This  is  the  proper  description  of  the  plasma 
in  a  thermionic  diode,  although  the  temperature  and  density  vary  across  the 
interelectrode  space.  The  plasma  is  optically  thin  to  radiative  decay  from 
upper  levels,  except  for  resonance  lines  (transitions  ending  on  the  ground 
state).  The  energy  lost  from  the  plasma  by  radiative  decay  can  then  be  pre¬ 
dicted  by  summing  the  energy  emitted  in  transitions  from  each  level  (omitting 
transitions  to  the  ground  state)  and  from  the.  continuum  (radiative  recombina¬ 
tion).  The  energy  lost  in  the  resonance  line  is  more  complicated  and  must 
be  considered  separately  for  each  condition  of  optical  thickness. 

The  radiated  energy  in  a  steady  state  plasma  is  provided  ultimately 
by  the  electrons.  This  is  because  in  the  model  that  has  been  adop;ed  the 
electrons  produce  excitation  and  ionization  of  the  atoms.  Hence,  the  radiative- 
energy  loss  rate  dW^/dt  is  the  electron  kinetic-energy  loss  rate.  This 
energy  must  be  balanced  by  the  discharge  (that  is,  the  cathode)  in  order  to 
maintain  the  plasma  in  steady  state. 

The  radiated  energy  can  be  simply  written  as 

dW  „  « 

dT  =  -  S  hvixi  N(<1)A(<I.P)  -  Ne  S  |  hv  3(p,v)dv  . 
p,q  p  o  1 


(l) 


The  first  term  on  the  right-hand  side  is  the  line  emission  resulting  from 
transitions  from  upper  level  q  to  lower  level  p  at  frequency  v  .  The 

rM 

density  of  atoms  in  level  q  is  N  q) ,  and  ACq^  is  the  Einstein  spontaneous 

3 

emission  rate.  The  second  term  is  continuum  emission  by  Ng  electrons/cm 
recombining  with  an  equal  density  of  ions.  The  velocity  dependent  recombina¬ 
tion  coeffidient  is  S(ptv)  arid  vp  is  the  frequency  of  the  transitions.  The 
frequencies  are  given  by 


and 


(2) 


where  Ep  and  £q  are  the  ionization  potentials  of  the  levels  p  and  q  . 
The  summations  are  over  all  levels  except  the  ground  state. 

The  rates  A  and  3  have  been  discussed  in  previous  work,  and  the 

populations  N(q)  are  given  by  the  computer  code  that  solves  the  rate  equa- 
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tions  for  the  atomic  levels.  '  ’  In  terms  of  our  earlier  notation, 


and 


N(p)  =  RQ(p)  +  R1Cp)NCl) 


ClNe  '  C2N(1)  ♦ 


(3) 

(4) 


where  represents  the  contributions  from  the  continuum  and  C£  describes 
the  energy  loss  resulting  from  excitation  of  the  ground  state  (p  =  1). 


The  computer  model  treats  53  levels  of  cesium,  of  which  the  first 
26  are  out  of  equilibrium  and  have  populations  given  by  Eq.  (3).  The  next  17 
levels  are  in  equilibrium  and  have  populations  given  by  the  Saha  equation. 
Still  higher  levels  have  negligible  populations  at  the  temperatures  considered 
here  and  are  neglected.  In  terms  of  these  levels  the  coefficients  and 
C2  are 


-26  53 

L  (Ep-  Eq)Ro(q)A(q,p)  +  £  (Ep-  Eq) Ng(q) A(q,p) 
q=3  CF=27 


+ 


J  (E  +  J£mv^)S(PtV)dv 


(5) 
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(6) 


26  26 

c2  =  I  £  <EP  '  Eq)Ri(£l)«<l.P> 
p=2  q=3 

These  coefficients  are  dependent  on  the  electron  density  and  temperature.  They 

have  been  calculated  from  solutions  of  the  rate  equations  using  the  code  pre- 
1  2 

pared  earlier  '  .  Complete  trapping  of  resonance  radiation  has  been  assumed 
in  determining  populations.  The  coefficients  C^  and  C2  are  like  the  ion¬ 
ization  rate  S(N  ,T  )  and  the  recombination  rate  cvCN  ,T  )  in  that  they  are 
e  e  e  e 

effective  rates  describing  the  result  of  many  transitions  between  levels. 

In  order  to  facilitate  practical  use  of  the  results,  the  following 
empirical  equations  were  developed  for  the  coefficients  a,  S,  and  C 2  . 


av.  =  (a  /N  )exp[E  / (kT  cm^/sec  (7) 

O  0  0  3  G 

S  =  SQ  exp[-Eb/kTe]  cm'Vsec  (8) 

C,  =  C,  exp[E  / (kT  )^]  watt-cm^  (9) 

1  1  r  c  e  J 

C2  =  Cl,  exp[-Ed/kTe]  watt  (10) 


The  empirical  constants  cf  ,  S  ,  C.1 ,  Cl,  E  ,  E,  ,  E  a:;d  E  are  functions  of 

electron  density  only  and  are  tabulated  in  Table  1.  All  the  temperature 
dependence  is  in  the  exponents.  The  multiplicative  constants  ,  So, 

C’  and  C2  are  in  the  units  given  in  Eqs.  (7)  -  (10),  while  the  exponen¬ 
tial  constants  E  ,  E,  ,  E  and  E,  are  such  that  kT  must  be  in  eV. 
abed  e 

Linear  interpolation  of  the  values  in  Table  1,  as  a  function  of  log  Ng  , 
will  reproduce  the  original  results  t-'  within  3%  in  the  range  1800  K  <  Tg 
<  3000°K. 

It  should  be  emphasized  again  that  resonance  radiation  is  not 
included  in  the  emitted  energy  given  by  C^  and  C2  . 

2.  Energy  Loss  Rates  in  a  Thermionic  Converter 

The  energy  loss  rates  described  in  the  previous  section  have  been 

4 

used  with  the  temperature  and  density  profiles  given  by  Reichelt  to  deter- 
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TABLE  I 


Coefficients  tor  use  in  Eqs.  (7)  -  (10)  for  the  calculation  of  effective 
coefficients  for  recombination,  ionization  and  radiative  energy  loss. 
Complete  trapping  (g  =  0)  of  resonance  radiation  has  been  assumed. 


|(Ne,cm"3) 

1A15 
a  x  10 

0 

E 

a 

S  x  105 

0 

Eb 

/  1A33 

C^  x  10 

E 

c 

C2  x  10 13 

Ed 

12.0 

26.00 

2.804 

1.378 

3.756 

69.03 

2.100 

3.665 

1.968 

12.1 

24.73 

2.880 

1.580 

3.735 

58.20 

2.200 

4.031 

1.966 

12.2 

23.32 

2.961 

1.783 

3.713 

48.35 

2.306 

4.425 

1.965 

12.3 

21.85 

3.048 

1.987 

3.688 

40.17 

2.409 

4.870 

1.967 

12.4 

20.40 

3.138 

2.191 

3.660 

33.35 

2.516 

5.378 

1.971 

12.5 

18.78 

3.235 

2.363 

3.630 

27.65 

2.624 

5.973 

1.976 

12.6 

17.19 

3.338 

2.494 

3.598 

22.95 

2.737 

6.640 

1.984 

12.7 

15.52 

3.447 

2.589 

3.564 

18.85 

2.852 

7.372 

1.992 

12.8 

13.87 

3.562 

2.630 

3.529 

15.32 

2.975 

8.180 

2.000 

12.9 

12.33 

3.681 

2.626 

3.492 

12.31 

3.100 

9.013 

2.008 

13.0 

10.95 

3.  MO 

2.576 

3.455 

9.829 

3.231 

9.860 

2.016 

13.1 

9.800 

3.917 

2.493 

3.419 

7.800 

3.361 

10.67 

2.021 

13.2 

8.868 

4.033 

2.388 

3.383 

6.230 

3.487 

11.42 

2.025 

13.3 

8.146 

4.141 

2.270 

3.348 

5.029 

3.608 

12.08 

2.028 

13.4 

7.675 

4.239 

2.149 

3.315 

4.126 

3.720 

12.63 

2.029 

13.5 

7.390 

4.331 

2.168 

3.285 

3.400 

3.825 

13.08 

2.028 

13.6 

7.369 

4.408 

1.925 

3.259 

2.865 

3.920 

13.39 

2.027 

13.7 

7.522 

4.479 

1.826 

3.235 

2.455 

4.001 

13.66 

2.025 

13.8 

7.933 

4.539 

1.742 

3.214 

2.128 

4.074 

13.83 

2.023 

13.9 

6.619 

4.589 

1.670 

3.196 

1.906 

4.135 

13.95 

2.020 

14.0 

9.484 

4.631 

1.609 

3.180 

1.738 

4.188 

14.04 

2.018 

14.1 

10.68 

4.665 

1.558 

3.168 

1.,  603 

4.230 

14.09 

2.016 

14.2 

12.30 

4.694 

1.517 

3.157 

1.499 

4.266 

14.12 

2.014 

14.3 

14.52 

4.718 

1.485 

3.148 

1.415 

4.297 

14.15 

2.013 

14.4 

17.20 

4.737 

1.457 

3.140 

1.350 

4.320 

14.17 

2.012 

14.5 

20.58 

4.751 

1.435 

3,135 

1.302 

4.341 

14.18 

2.011 

14.6 

24.71 

4.763 

1.418 

3.130 

1.268 

4.357 

14.19 

2.010 

14.7 

30.52 

4.773 

1.405 

3,127 

1.238 

4.369 

14.20 

2.009 

14.8 

37.55 

4.780 

1.393 

3.124 

1.213 

4.379 

14.20 

2.008 

14.9 

46.41 

4.786 

1.384 

3.121 

1.194 

4.387 

14.21 

2.008 

15.0 

57.37 

4.792 

1.377 

3.119 

1.179 

4.395 

14.21 

2.007 
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mine  the  actual  energy  lost  by  radiation  in  a  typical  case.  The  energy 

emitted  in  the  resonance  lines  has  also  been  included  by  using  an  approximate 

solution  for  the  radiation  transport  equation.  The  lines  are  optically  thick 

with  considerable  energy  emission  in  the  wings.  The  total  emission  has  been 

5 

obtained  using  Holstein's  escape  probabilities,  with  pressure  broadened 
line  shapes  in  the  wings. 

Details  of  the  calculation  and  the  results  have  been  presented  at 
the  Physical  Electronics  Conference,  Cambridge,  Massachusetts,  1967.  The 
paper  is  included  in  the  Proceedings  and  is  reproduced  as  Appendix  1  of  this 
report. 


2.3  Plasma  Transport  Equation 

The  differential  equations  governing  the  transport  of  particles, 

momentum,  and  energy  in  the  plasma  between  the  electrodes  of  a  thermionic 

converter  have  been  formulated  in  a  more  complete  fashion  than  has  been  done 

previously.  Two  specific  advances  over  earlier  work  are  the  inclusion  of 

£ 

more  exact  coefficients  for  electron  transport  and  for  ion  production  and 
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radiative  energy  loss  in  the  plasma  volume. 

In  this  section  these  equations  are  formulated  and  cast  in  a  form 

suitable  for  computer  integration. 

—2  —1  6 
The  electron  flux,  cm  -sec  ,  is  given  by  the  equation 


r  =  -  u 


kT  dN  N  dT 

NeE  +  -r  ~i  +  T^-v2) 


-3  -1 

The  electron  kinetic  energy  flux,  erg-cm  -sec  ,  is 


N  dT 

°  =  5  kTe  r  -  f  kTe  -  5‘ }  k  if 


(11) 


(12) 


where  Ng  is  the  electron  (ion)  density,  T  is  the  electron  temperature, 
k  is  Boltzmann's  constant,  E  is  the  electric  field,  and  is  the 
electron  mobility.  These  are  standard  equations,  including  the  effect  of 
both  electron  pressure  and  thermal  diffusion. 


*  The  parameters  pe  ,  §  and  Q  are  functions  of  Tg  and  the  degree  of 
ionization. 
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The  equation  for  ion  flux,  cnf^-sec"1.  corresponding  to  Eq.  (11)  is^ 

N_  dT 


r  =  p 

p  p 


N  E 
e 


kT  dN  N  dT  N 

-£  _£  _  _£  k  -H  _  kT  -S  k 
e  dx  e  dx  i  e 


5,1 

dxj 


(13) 


T 

where  T  is  the  ion  temperature,  p  is  the  ion  mobility,  and  K.  is  the 
P  P  1 

ion  thermal-diffusion  constant.  Since  thermionic  converter  plasmas  are 

characterized  by  relatively  low  degrees  of  ionization,  the  term  in  Ref.  (8) 

which  accounts  for  transfer  of  directed  electron  momentum  to  the  ions  has  not 

been  included. 


In  addition  to  Eqs.  (11)  -  (13)  continuity  equations  for  the  para¬ 
meters  T  ,  Q  and  r  are  required.  We  will  assume  that  the  electron  flux  is 
constant  and  given  by  the  measured  current  so  that 


(14) 


The  electron  kinetic-energy  flux  is  decreased  by  the  energy  costs  of  ion 
production  and  radiative  loss,  and  increased  by  Joule  heating,  expressed  by 


^-E 

dx 


-  R  -  TeE 


(15) 


where  1^  is  the  first  ionization  energy  of  cesium  and  -R  is  the  radiative 
energy  loss.  The  ion  production  rate  is 


=  5  Vl  -  (16) 

where  S  and  are  the  effective  three-body  ionization  and  recombination 
coefficients  (functions  of  N  and  T  ) ,  and  Nj  is  the  density  of  ground- 
state  cesium  atoms. 

Equations  (11)  -  (16)  are  the  primary  equations  in  this  formulation. 
If  we  make  the  assumptions  that  the  heavy  particle  temperature  varies  linearly 
between  the  electrode  temperatures,  and  that  the  ideal  gas  law  is  valid  in 
the  plasma,  we  can  also  write  the  equations 

dT 

— J=-(TE-TC)/L  (17) 
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Na=kf 


where  Tg  and  Tg  are  the  emitter  and  collector  temperatures,  respectively, 

P  is  the  cesium  vapor  pressure,  and  L  is  the  electrode  spacing. 

Returning  to  Eq.  (12)  the  parameters  ,  £  and  £  are  derived 
in  Ref.  (6)  under  the  assumption  that  electron-electron  effects  are  negligible. 
This  assumption  is  valid  for  low  degrees  of  ionization.  In  Ref.  (6)  the  para¬ 
meters  §  ,  C  «  and  the  effective  collision  frequency  v  ^  (including  both 
electron-atom  and  electron-ion  effects)  are  presented  graphically  as  functions 
of  T  and  the  degree  of  ionization  N  /N  .  The  mobility  u  is  related  to 

C  c  d  v 


^e  mv 


where  m  is  the  electron  mass.  The  graphical  data  of  Ref.  (6)  have  been 
converted  to  Tables  2,  3  and  4  for  our  use. 

If  we  make  the  further  assumption,  based  on  a  low  degree  of  ion¬ 
ization,  that  the  ion-ion  collision  frequency  is  much  less  than  the  ion-atom 

collision  frequency,  the  ion  thermal-diffusion  constant7  can  be  taken  to  be 
T 

K.  -  0.18.  The  ion  mobility  is  related  to  the  value  at  standard  conditions 

by9  N  /  T  f 

V-  „  =  ^ktU2)  (20) 


where  N  =  2.69  x  IQ19  atoms-cm-3,  T  =  273°K  and  \i  =  0.12  cm2  volt-1-sec-1.3,9 
0  T  0  0 

The  radiative  energy  loss  term  in  Eq.  (15)  can  be  written^ 

-  R  =  C.N2  +  C9N.  +  5.16  x  10~14  G  N. e-1 '432^kTe  (21) 

1  e  i  i  1 

where  G  is  the  geometrical  factor  for  loss  of  resonance  radiation.  In  terms 
of  the  average  electrode  emissivity  e  ,  plasma  length  L  ,  and  radius  R 

G  =  — +  — •  .  (22) 


The  coefficients  and  C£  are  functions  of  Ng  and  Tg  and  are  given 
by  Eqs.  v.9)  and  (10).  The  coefficients  S  and  are  given  by  Eqs.  (7) 
and  (8) . 
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TABLE  II  -  v  /N 

,  (cm^  sec 

1  > 

loo 

Ne 

V°K 

Ji0 

N 

a 

1750 

2000 

2250 

'2500 

2750 

3000 

-3.000 

4.64(- 

7) 

4 .72( — 7) 

4.87 ( —7) 

5 . 19  ( —7) 

5 . 53( —7) 

5.90C-7) 

-2.699 

6.09(- 

7) 

6 . 19 ( —7) 

6.39C-7) 

6.71C-7) 

7.03C-7) 

7.50(-7) 

-2*301 

1.000( 

-6) 

1.016C-6) 

1.049 (-6) 

1.083C-6) 

1.118C-6) 

1.173(-6) 

-2.000 

v — / 

CM 

vO 

• 

-6) 

1.615C-6) 

1.615C-6) 

1.642C-6) 

1.695(-6) 

1.750(-6) 

-1.699 

2.74(- 

6) 

2.65C-6) 

2.6K-6) 

2.63C-6) 

2.66(-6) 

2.70C-6) 

-1.301 

5.62(- 

6) 

5 . 27 ( —6) 

5.03C-6) 

4.87(-6) 

4.72(-6) 

4.64C-6) 

-1.000 

1.033C 

-5) 

9.08C-6) 

8.39C-6) 

7 .87( —6) 

7.38C-6) 

7 .03( —6) 

TABLE  III  -  C  , 

Dimensionless 

loo.  „ 

N 

e 

Te.°K 

J10 

N 

a 

1750 

2000 

2250 

2500 

2750 

3000 

-3.000 

6.71 

5.48 

4.71 

4.20 

3.85 

3.73 

-2.699 

7.01 

5.82 

5.01 

4.58 

4.24 

4.11 

-2.301 

7.57 

6.46 

5.73 

5.30 

5.09 

5.13 

-2.000 

8.29 

2.27 

6.63 

6.33 

6.20 

6.37 

-1.699 

9.32 

8.42 

7.91 

7.74 

7.78 

8.08 

-1.301 

11.24 

10.55 

10.26 

10.34 

10.60 

11.07 

-1.000 

13.03 

12.56 

12.43 

12.65 

12.99 

13.50 

TABLE  IV  ..  g  , 

Dimensionless 

.  \ 

Ne 

Te,°K 

lo>i0 

N 

a 

1750 

2000 

2250 

2500 

2750 

3000- 

-3.000 

2.46 

2,23 

2.04 

1.89 

1.78 

1.70 

-2.699 

2.50 

2.27 

2.09 

1.94 

1.84 

1.77 

-2.301 

2.58 

2.36 

2.19 

2.07 

1.98 

1.93 

-2.000 

2.68 

2.48 

2.32 

2.22 

2.15 

2.12 

-1.699 

2.81 

2.64 

2.51 

2.43 

2.39 

2.38 

-1.301 

3.05 

2.91 

2.82 

2.78 

2.77 

2.80 

-1.000 

3.25 

3.16 

3.09 

3.08 

3.09 

3.13 

-  9  - 


2 

The  pressure  P  ,  dyne-cm  ,  in  Eq.  (18)  is  related  to  the  cesium 
reservoir  temperature,  T^  ,  by  the  Langmuir  equation 

P  =  1.332  x  103  Tr  _1,35  1011-0531-4041*°/TCs 
Ls 

The  ground  state  density  and  total  atom  density  are  related  by  =  IlN^  , 
where  n  is  the  partition  function  for  cesium 


i=l 


I./kT 
l  e 


(24) 


The  statistical  weights  and  ionization  energies  of  the  bound  levels  of  cesium 
are  g.  and  I.  ,  respectively. 

It  is  more  convenient  to  use  the  Einstein  relation 


u-  - 


eD 

kT 


(25) 


in  rewriting  Eqs.  (1D-(13).  After  some  algebra,  the  transport  equations  can 
be  transformed  to  the  form 


dTe  5kTer  -  Q 


dx  (r  -  52)kN  D 

Pi 


e  e 


(26) 


dNf 

dx 


T  +  T 
e  p 


/  r  t  r  n 

!  T  ^ 

i  dT 

dT  ' 

(—  +  _E  _fi  +  _e 
\d  t  d  t 

1  +  kt 

to12  +  (5  - 

3/2) 

1  e  e  p  e 

i  1 

1 

_ 

(27) 


dr  , 

T2  --  SN  N,  -  <*,pr 
dx  el  3  e 


dV  = 
dx 


=  -  E  =  - 


kTp 

eN  D 
e  p 

dfl- 


dN  D  dT 

r  +  d  ~  n  — ^ 

p  p  dx  T  e  dx 
.  P 

dr 


.  =  -  I  t-2  -  R  -  TEe 
dx  p  dx 


(28) 


(29) 


(30) 


dT 

_e  _ 

dx  ~ 


(TE-  tc)/l 


(31) 
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3.0  DIODE  EXPERIMENTS 

3.1  Objectives 

The  objectives  of  the  spectroscopic  work  are  to  provide  further 
data  on  the  spatial  variation  of  electron  temperature  and  density  across  the 
interelectrode  gap  of  a  converter  for  experimental  verification  of  the  non¬ 
equilibrium  calculations  and  to  explore  features  of  the  spectrum  that  may  be 
used  for  refined  diagnostic  measurements.  In  addition,  ionization  cross 
sections  have  been  measured,  with  particular  emphasis  on  the  threshold  behavior. 


3.2  Principles  of  Spectroscopic  Measurements 

3.2.1  Recombination  Intensities.  Electron  temperatures  have  been  deter¬ 
mined  from  measurement  of  the  continuum  emitted  by  radiative  recombination  to 
the  6P  states  of  cesium.  The  reaction  is 

Cs+  +  e  -  Cs(6P)  +  hv  (32) 

The  6P  is  commonly  used  because  it  is  usually  the  strongest  continuum  and  it 
is  in  a  region  of  the  spectrum  that  is  relatively  free  from  lines.  The  inten¬ 
sity  distribution  of  the  continuum  enables  the  determination  of  the  electron 
temperature  using  only  relative  intensity  measurements.11 

o 

The  intensity  distribution  has  its  threshold  at  about  5080  A  and 
is  given  by 

I(X)dX  =  hvN  N.v  a  (v)f(v)dv  (33) 

3 

where  I(X)dX  is  the  emission  rate  (erg/cm  -sec)  in  the  wavelength  interval 

dX  about  X  ,  corresponding  to  electrons  in  the  velocity  interval  dv  about  v  , 

a  is  the  velocity-dependent  cross  section  for  recombination  to  the  6P  level 

and  f(v)dv  is  the  fraction  of  electrons  in  dv  about  v  .  The  most  recent 

12 

values  for  the  cross  section  are  due  to  Agnew  and  Summers 


cWv)  = 


3.74  x  10 


-6 


cm 


(34) 


These  compare  with  earlier  measurements  by  Mohler 
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of 
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(35) 


where  X  is  the  wavelength  in  Angstroms.  The  electron  velocity  v  is  in 
cm/sec  in  both  equations. 


The  Agnew  and  Summers  cross  section  is  more  correct,  since  it  was 
obtained  from  measurements  of  an  equilibrium  system  (isothermal).  Their 
cross  section  and  a  Maxwell  distribution  of  velocities  in  Eq.  (33)  gives 


I(X)  =  3.74  x  10"6  N  N.  2m 

e  i  tt 


-E/kT 


or 

4n[x3I(\)]  =  K  -  3/2  4nT  -  E/kT  , 


where  K  is  a  constant  and  E  is  the  energy  above  the  series  limit.  Since 
the  temperature  is  constant  during  recording  of  the  spectrum,  a  plot  of 
Jtntx' I(X)  ]  vs  E  gives  a  straight  line  with  a  slope  of  -1/kT.  If  the 
Mohler  cross  r.ection  is  used,  the  logarithmic  equation  is 


4n[Xl(X)]  =  M  -  3/2  AnT  -  E/kT  t  (37) 

where  M  is  a  different  constant.  The  only  important  change  here  is  the 

power  of  X  multiplying  I(x)  on  the  left  side.  Over  the  small  wavelength 

o 

range  observed  in  typical  experiments  (4000  -  5000  A)  the  Mohler  cross  section 
gives  temperatures  greater  by  about  15%.  Equation  (36)  is  preferred. 

The  recombination  continuum  can  also  be  used  to  determine  electron 
density.  In  this  case  absolute  measurement  of  the  intensity  I(X)  must  be 
made.  Equation  (33)  is  then  solved  for  using  the  temperature  from 

the  relative  shape  of  the  continuum  to  obtain  F(v).  The  ion  and  electron 
density  can  be  assumed  equal  in  most  converters,  so  Ng  is  determined. 

Mohler* s  cross  section  in  this  case  gives  densities  about  35%  higher  than 
those  obtained  from  Agneiv  and  Sommer's  cross  section.  The  absolute  measure¬ 
ment  is  rather  difficult  to  perform  and  is  subject  to  absorption  errors. 

As  a  result,  the  more  easily  performed  and  less  sensitive  line  shape  measure¬ 
ment  is  considered  more  reliable. 


-  12  - 


/ 


3.2.2  Line  Shapes.  Line  shapes  give  densities  when  the  lines  are  pre¬ 
dominantly  Stark  broadened  as  opposed  to  other  broadening  mechanisms.  In 
such  a  case  the  width  is  nearly  directly  proportional  to  the  electron  density. 
The  widths  have  been  worked  out  in  detail^' ^  and  verified  experimentally.^ 
The  basic  mechanism  that  results  in  broadening  is  the  shift  of  the  levels  of 
the  radiating  atom  due  to  the  electric  field  (Stark  effect)  of  the  electrons 
and  ions  of  the  plasma.  The  shift  is  averaged  over  the  distribution  of  elec¬ 
tron  and  ions  around  the  radiating  atom  and  results  in  a  broadened  (and 
shifted)  line.  The  theory  has  been  developed  for  isolated  lines  so  that  the 
broadening  must  not  be  so  large  that  lines  overlap. 

In  practice  the  fundamental  series  lines  are  most  suitable  for  Stark 

13  15  3 

broadening  measurements.  Over  a  density  range  from  10  -  10  electrons/cm0 

the  lines  are  isolated  and  sufficiently  wide  that  their  shape  can  be  deter¬ 
mined  with  low  resolution  spectrometers.  The  shapes  are  only  slightly  tem¬ 
perature  dependent  (a  few  percent  for  cesium  over  temperatures  from  2000°K  to 
5000°K)  and  so  the  temperatures  need  not  be  known  to  determine  electron  den¬ 
sity.  The  density  is  determined  in  practice  by  measuring  the  widths  of 
several  fundamental  series  lines  at  half  maximum  intensity  and  comparing 
with  calculated  values. 

3.3  Experimental  Procedure 

The  diode  was  run  at  a  constant  body  and  emitter  temperature  and 
with  constant  gas  flow  to  the  collector  seal.  The  cesium  density  was  adjusted 
by  varying  the  reservoir  temperature.  This  was  done  from  low  to  high,  always 
keeping  it  enough  lower  than  all  other  sections  so  that  it  was  the  controlling 
element.  A  factor  of  two  increase  in  pressure,  which  required  a  rise  of  about 
25°C,  would  stabilize  in  15  minutes.  The  temperature  was  controlled  by  vary¬ 
ing  the  current  through  the  diode.  In  this  case  stabilization  was  almost 
instant  and  a  run  could  be  repeated  after  five  minutes.  The  diode  current 
control  circuit  is  shown  in  Fig.  1. 

Stability  was  verified  by  peaking  the  monochrometer  manually  on  any 
of  the  lines  of  interest  and  running  for  five  minutes  with  zero  drift. 

Emitter  temperature  was  monitored  continuously  with  a  two  color  optical  pyro¬ 
meter.  This  device  takes  the  ratio  of  the  color  temperature  at  two  different 
bands  of  wavelengths  in  the  visual  spectrum.  This  helps  to  correct  for 
variations  in  emissivity  of  the  emitter.  The  area  of  the  emitter  viewed  was 
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1  Optical  bench  layout 


a  spot  0.1“  in  diameter  on  the  center  line  of  the  sleeve  immediately  behind 
the  emitter  button.  The  weld  line  where  the  button  was  attached  to  the  sleeve 
was  used  as  the  reference  for  one  edge  of  the  area  viewed.  This  gave  con¬ 
tinuous  monitoring  without  any  influence  from  the  plasma  except  at  the  lowest 
densities  at  which  measurements  could  be  made. 

Calibration  for  the  measurement  of  relative  and  absolute  line 
intensities  was  made  by  replacing  the  diode  with  a  tungsten  strip  lamp  cal¬ 
ibrated  as  a  standard  of  spectral  radiance  by  NBS  (Fig.  2).  This  substitution 
procedure  obviated  the  need  for  corrections  dm.  to  reflection  losses,  variation 
in  solid  angles,  etc.  All  line  intensity  measurements  and  the  majority  of 
the  line  broadening  measurements  were  made  using  a  0.5  meter  scanning  mono¬ 
chrometer  (Jarrell-Ash  Model  82-000)  equipped  with  a  15,000  lpi  grating.  For 
the  measurements  of  forbidden  line  intensities  and  shapes  and  resonance  line 
shapes  a  1.0  meter  scanning  monochromator  (Jarrell-Ash  Model  178-420)  was 

employed.  This  instrument  was  equipped  with  a  30,000  lpi  grating  blazed  at 

o 

lp.  A  bandpass  of  better  than  0.03  A  (3rd  order)  was  demonstrated  with  this 
instrument.  This  monochromator  was  also  used  to  check  the  line  shapes  measured 
with  the  lower  resolution  spectrometer.  Little  or  no  correction  for  instru¬ 
ment  width  was  required. 

Since  one  of  the  objectives  of  this  study  was  to  measure  recombina- 
tion  temperatures  at  low  electron  densities  (<  10  /'cm  ),  considerable 
effort  was  expended  to  maximize  the  signal-to-noise  ratio  of  the  system. 

Phase  sensitive  detection  was  tried  but  the  high  time  constants  which  are 
required  in  this  technique  made  it  unsuitable  for  our  purposes.  The  most 
successful  technique  was  the  use  of  a  magnetic  lens  (Fig.  3).  This  consists 
simply  of  a  200  ampere-turn  coil  coaxial  with  the  photocathode  of  the  EMI 
9558  A  photomultiplier  tube.  This  had  the  effect  of  reducing  the  system  noise 
by  a  factor  of  100  while  increasing  the  signal  about  15%.  This  arrangement 
provided  the  requisite  sensitivity  for  the  measurement  of  the  low  light  levels 
encountered. 

3.4  Temperature  and  Density  Spatial  Distributions 

The  external  optical  system  allowed  the  interelectrode  plasma  to 
be  focused  on  the  entrance  slit  of  the  monochromator  with  a  magnification 
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FIG.  2  Photomultiplier  detection  system 
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of  5.5X.  The  monochromator  itself  was  mounted  on  tracks  which  allowed  it  to 

be  moved  transverse  to  the  optical  axis  with  a  precision  of  0.1  mm.  This 

allowed  spatially  resolved  spectral  measurements  to  be  made.  In  this  manner 

the  temperature  and  electron  density  profiles  shown  in  Figs.  4  and  5  were 

obtained.  As  in  all  of  this  work,  the  densities  were  measured  by  observing 

o 

the  half  width  of  the  8F  -  5D  transition  at  6628  A.  This  half  width  was 

related  to  the  density  in  the  manner  described  above.  Temperatures  were 

measured  by  measuring  the  intensity  variation  of  the  recombination  continuum 

o  oo 

as  described  above.  Four  points  at  100  A  intervals  from  5000  A  to  4700  A 

were  measured  for  the  data  reported  here.  Later  measurements  increased  the 

o 

lower  wavelength  range  to  4300  A  but  did  not  alter  the  conclusion  significantly. 
(These  data  are  too  incomplete  for  inclusion  here.)  The  recombination  data 
was  reduced  by  computer.  At  the  same  time  the  temperatures  were  obtained, 
density  data  were  also  generated  because  an  absolute  calibration  of  the 
optical  system  was  available.  These  density  measurements  were  not  in  satis¬ 
factory  agreement  with  the  densities  measured  using  line  width.  The  latter 
densities  are  considered  more  accurate,  since  less  precise  measurements  are 
required  in  obtaining  them  and  the  line  broadening  technique  is  well  established. 
The  difference  appears  to  be  due  to  extraneous  radiation  in  the  gap,  but  time 
did  not  permit  a  definitive  study  of  this  problem. 

3.5  Resonance  Line  Shapes 

Resonance  lines  have  been  observed  in  emission  under  a  variety  of 
plasma  conditions.  They  are  self-absorbed  because  of  the  cold  gas  outside  the 
plasma  region.  However,  they  have  a  characteristic  asymmetry,  as  shown  in 
Fig.  6.  There  appears  to  be  a  shift  of  the  emission  line  relative  to  the 
absorption  line  of  the  cold  gas  and  also  an  asymmetry  of  either  the  absorption 
or  emission  line,  or  both. 

The  resonance  absorption  lines  have  also  been  measured  by  shining 
a  standard  lamp  through  the  cold  cesium  gas  with  the  diode  turned  off. 

Figures  7  and  8  show  the  shapes  of  the  6S  -  7P  doublet  lines  as  pressure 
is  varied.  The  lines  have  very  surprising  asymmetry.  The  cause  is  not  clear, 
but  is  probably  due  to  small  collisional  broadening  of  the  line.  The  effect 
is  magnified  in  an  optically  thick  gas  as  the  two  wings  of  a  line  have  dif¬ 
ferent  absorption  and  saturate  at  different  optical  thicknesses. 
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FIG.  6  Shape  of  the  Cs-4593  A  line  in  emission. 
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3.6  Forbidden  Line  Intensities 

On  the  low  wavelength  side  of  the  Stark  broadened  fundamental  series 

lines  (5D  -  nF) ,  relatively  weak  "forbidden"  lines  of  the  6F  -  5D  series  are 

observed.  The  forbidden  and  allowed  lines  are  well  separated  for  transitions 

to  the  lower  F  levels  and  become  stronger  and  closer  together  for  trans- 

o 

itions  to  higher  F  levels.  The  6G  -  5D_  transition  (7269  A)  at  an  electron 
Id  -3 

density  of  1.5  x  10  cm  is  shown  in  Fig.  9.  The  allowed  6F  -  5D line 

at  7279  A  has  a  peak  intensity  well  above  the  top  of  the  paper.  The  same  line 

at  two  higher  electron  densities  is  shown  in  Fig.  10.  The  electron  densities 

o 

have  been  determined  from  Stark  broadening  measurements  of  the  6628  A  (8F  - 
line. 

The  forbidden  line  has  definite  structure  in  the  center  and  very 
sharp  fall-off  in  the  wings.  The  center  appears  to  be  made  up  of  two  lines 
slightly  displaced,  giving  a  double  peak*  As  the  density  is  increased,  the 
line  becomes  broader,  the  total  intensity  increases,  and  the  double-peak  struc¬ 
ture  is  lost. 

Considerable  data  were  taken  of  the  total  intensity  of  the  forbidden 

o 

line  as  the  density  was  varied.  These  results  for  the  7270  A  (6G  -  5Q5/2) 

lines  are  shown  in  Fig.  11.  The  intensity  of  this  line  divided  by  the  inten- 

0 

sity  of  the  allowed  6034.09  A  (lOSj/2  ~  ^p3/2^  line  is  plotted  against  electron 

density.  Division  by  the  intensity  of  the  6034.09  A  line  has  been  done  to 

eliminate  temperature  effects  that  change  intensities  as  density  is  varied. 

0 

The  upper  state  of  the  6034.09  A  line  (lOS^l  lies  very  near  the  upper  state 
of  the  forbidden  line  (6G) .  Therefore,  population  changes  as  the  temperature 
is  varied  are  the  same  for  both  lines  and  the  ratio  should  be  independent  of 
temperature. 

Figure  11  shows  a  linear  increase  of  total  intensity  with  density 
until  about  4.2  x  10^  cm  ^  .  Thereafter,  the  intensity  continues  to  increase 
linearly  but  at  a  lower  rate.  The  density  at  this  rather  abrupt  change  in 
slope  is  approximately  where  the  double-peak  nature  of  the  forbidden  line 
disappears . 

The  forbidden  lines  are  present  because  of  the  electric  field 
(Stark  effect)  of  the  ions  and  electrons  in  the  plasma.  These  fields  destroy 
atomic  degeneracy  of  levels  and  shift  the  energies  of  the  levels.  The  theory 
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FIG.  9  Forbidden  line  profile  of  7269  A  (6G-5Dg/2)  on  wing  of  the 
allowed  7279  A  (6F  -5D5/2)  as  measured  with  the  high  resolu¬ 
tion  instrument.  N  =  1.5  x  10^  cnT^;  T  »  2000°K. 

e  e 
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A  A  (A) 

FIG.  10  The  7269  A  forbidden  profile  showing  the  effect  of  electron 

density.  Ne  =  4.85  x  1014  for  solid  line  and  2.20  x  1014  cnr^ 

for  broken  line.  T  «  2000°K. 
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FIG.  11  Intensity  ratio  of  Cs  7269  A  (6G  -  5D5/2)/Cs  6034  A  (lOS^  -  6P3/2) 
vs  Ne  measured  by  broadening  of  Cs  6628  A  C8F  -5D5/2). 
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is  quite  analogous  to  the  Stark  broadening  theory  of  lines,  and  the  forbidden 
lines  are  a  second  manifestation  of  the  plasma  forces. 

Some  theoretical  work  has  been  done  on  forbidden  line  intensities 
in  plasmas,^  but  only  for  isolated  lines.  Isolated  lines  do  not  overlap 
other  lines;  in  an  element  like  cesium  this  means  that  the  upper  level  of  the 
line  is  not  degenerate  with  any  other  level.  Unfortunately,  the  transitions 
discussed  here  are  not  isolated  because  the  upper  level,  6G,  is  degenerate 
with  6H  and  61.  Therefore,  hydrogenic  theory  (levels  degenerate  in  i  )  must 
be  used  rather  than  isolated  line  theory.  This  has  not  yet  been  worked  out. 

The  double  peaks  of  the  5D  -  6G  lines  are  probably  a  result  of  the 
degeneracy  of  the  6G  and  6H  levels.  Splitting  into  double  peaks  is  character¬ 
istic  of  Stark  broadening  of  degenerate  levels  and  is  very  pronounced  in 
hydrogen.  The  disappearance  of  the  double  peaks  as  density  is  increased  is 
quite  likely  a  result  of  the  broadening  of  each  of  the  two  lines  result  in 
the  observation  of  one  peak. 

The  data  of  Figs.  9  -  11  is  the  first  quantitative  data  of  forbidden 
non-isolated  line  intensities  induced  by  plasma  effects.  Theoretical  explan¬ 
ation  of  the  intensities  and  shapes  will  be  a  critical  and  needed  test  of 
Stark  broadening  theory.  Moreover,  the  mere  presence  of  forbidden  lines  is 
a  rough  indication  of  density,  and  their  observation  can  serve  as  a  quick 
diagnostic  measurement.  Quantitative  density  determinations  can  be  obtained 
from  curves  such  as  Fig.  11,  and  such  measurements  serve  as  an  alternative 
to  Stark  broadening  measurements. 
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4.0  ELECTRON  IMPACT  IONIZATION  CROSS  SECTION  IN  CESIUM 

4.1  Present  Status 

The  interpretation  of  many  plasma  experiments  using  cesium  (for 

instance  thermionic  energy  converters,  ion  propulsion  devices,  and  vapor 

release  in  the  atmosphere),  rests  upon  a  full  understanding  of  the  various 

mechanisms  of  ionization.  The  objective  of  our  work  has  been  to  determine 

the  cross  section  for  ionization  of  cesium  atoms  by  electrons.  Optical 

absorption  measurements  by  Kratz1*  have  given  the  series  limit  in  cesium  at 

3.8926  eV.  Electron  beam  measurements  agree  with  this  value,  but  are,  of 

course,  not  as  accurate.  The  efficiency  of  ionization  (in  arbitrary  units) 

of  the  cesium  atom  by  electron  bombardment  was  first  measured  by  Tate  and 

Smith  °  with  energies  ranging  from  threshold  to  700  eV.  They  discriminated 

between  singly-  and  multiply-charged  ions  with  a  mass  spectrometer.  The 

absolute  magnitude  of  the  cross  section  was  first  determined  in  the  crossed- 

19 

beam  experiment  of  Brink,  followed  by  more  measurements  in  another  apparatus 

20 

of  the  same  type  by  McFarland  and  Kinney.  In  the  crossed-beam  experiments 

the  data  were  obtained  for  energies  between  30  and  500  eV  in  steps  of  100  eV 

above  100  eV.  These  high-energy  results  were  later  supplemented  by  the 

2 1 

threshold  and  low-energy  measurements  of  Heil  and  Scott  1  (up  to  50  eV)  and 

oo 

Korchevoi  and  Przonski  (up  to  25  eV) ,  who  all  used  systems  in  which  cesium 
vapor  was  in  thermodynamic  equilibrium  with  a  liquid-phase  reservoir. 

19  20  22  21 

The  ratio  between  the  highest  r  *  and  lowest  maximum  cross- 
section  values  reported  is  about  a  factor  of  two.  The  reasons  for  this 
discrepancy  are  not  completely  understood.  On  this  basis,  we  found  it  desir¬ 
able  to  make  a  new  determination  of  the  ionization  cross  section  from  threshold 
to  100  eV  in  order  to  get  an  overlap  between  the  low-  and  high-energy  measure¬ 
ments  discussed  above, 

18 

We  used  a  Tate  and  Smith  -type  apparatus  modified  to  include  the 

23 

feature  of  retarding  potential  difference  (RPD)  in  the  electron  gun.  It 
is  crucial  that  the  number  density  of  cesium  atoms  in  the  collision  j|gion 
be  accurately  known.  Measurements  with  a  surface  ionization  detector  yielded 
density  values  that  agreed  with  those  calculated  from  the  formulas  of 
Taylor  and  Langmuir  within  +  3%. 
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4.2  Experimental 

4.2.1  General  -  The  major  components  of  the  experiment  are  shown  in  Fig.  12. 
The  electron  gun  and  ionization  region  (to  be  described  later)  are  located 

in  the  middle  of  a  vertically  hanging  tube  with  a  pool  of  cesium  at  the  bottom. 
The  cross  in  the  wide  pump-out  tube  contains  an  ionization  gauge  for  measur¬ 
ing  background  pressure  during  baking  and  pumping  and  a  surface  ionization 
detector  for  determination  of  cesium  density.  The  vacuum  parts  are  all  made 
from  stainless  steel  and  copper.  During  bake-out,  the  vertical  tube,  valves 
V.  and  V0,  and  the  Vaclon  pump  to  the  left  are  all  heated  to  300°  -  400°C 
while  pumping  is  done  with  the  other  Vaclon  pump.  Cesium  °  is  then  evaporated 
into  the  pool  through  a  small  side-arm  (not  shown  in  the  illustration). 

Under  normal  operation  the  high-vacuum  valves  are  closed,  and  the  ‘'clean" 

Vaclon  pump  to  the  left  in  Fig.  12  acts  on  the  experimental  volume  through 
a  small-diameter  copper  tubing  from  the  cesium  pool.  Part  of  this  tube  is 

water-cooled  to  about  35°C  in  order  to  prevent  cesium  from  entering  the  pump. 

-9 

An  end  vacuum  of  10  torr  is  easily  obtained  after  the  bake-out  procedure. 

The  apparatus  shown  in  Fig.  12  can  be  moved  up  and  down  as  one  unit,  which 
facilitates  the  lowering  of  the  tube  into  a  solenoid.  Circulating  hot  air 
keeps  the  temperature  of  the  tube  and  the  valve  V^  at  about  12D°C,  thereby 
avoiding  undesirable  cold  spots  where  the  cesium  could  condense.  The  temper¬ 
ature  of  the  Cs  pool  is  kept  constant  to  within  +  0.05°C  by  means  of  a 
thermostatically  controled  closed  loop  of  oil.  Notice  that  Cs  atoms  effuse 
out  of  the  lower  orifice  of  the  small-diameter  pumping  tube  with  the  temper¬ 
ature  TCs  ,  since  this  part  of  the  tube  is  immersed  in  the  oil  bath.  The 
length  immersed  was  actually  helix-shaped,  thereby  making  the  constant- 
temperature  region  much  longer  than  indicated  in  Fig.  12.  It  is  also  impor¬ 
tant  that  the  diameter  of  the  pump-out  tube  be  much  smaller  than  the  diameter 
of  the  cecium  pool. 

4.2.2  Cross  Section  Apparatus  •  The  electron  gun  in  Fig.  13  was  constructed 

23 

according  to  the  RPD  principle  of  Fox  _et  al.  All  electrodes  were  made 
from  copper.  Electrons  were  emitted  from  an  indirectly  heated  Philips  cathode 
of  type  A.  The  emission  current  from  this  cathode  in  the  presence  of  cesium 
was  found  to  be  very  stable,  so  that  stabilization  of  the  current  by  elec¬ 
tronic  means  was  unnecessary. 
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Vaclon  pump 


FIG.  12  Schematic  diagram  of  vacuum  system  and  cesium  reservoir.  Surface 
ionization  detector  and  ionization  gauge  are  mounted  on  the  two 
flanges  above  the  ionization  region. 
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FIG.  13  Cross  section  apparatus.  The  spacing  between  electrodes  in  the  gun  was 
1  mm.  The  height  of  the  aperture  in  the  retarding  electrode  was  0.2  mm. 
The  length  of  the  ion  collector  was  4  =  26.6  iron  .  The  size  of  the  gaps 
between  this  collector  and  the  guard  plates  is  strongly  exaggerated. 
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The  potentials  indicated  below  the  gun  in  Fig.  13  are  all  with 
respect  to  the  cathode.  All  voltages  were  supplied  from  batteries.  The  narrow- 
slit  retarding  electrode  is  at  a  potential  -V^  .  Only  electrons  with  energies 
above  will  enter  the  acceleration  region  and  gain  an  energy  before 

entering  the  collision  region  between  the  parallel  plates.  A  small  change 
AVpj  in  the  retarding  potential  gives  a  change  Al  in  the  electron  current 
and  a  corresponding  change  AI+  in  the  ion  current.  If  only  singly  charged 
ions  are  produced,  the  cross  section  can  be  determined  from  the  relation 


a 


1 

nZ  5 


(38) 


where  n  is  the  cesium  atom  number  density  and  &  is  the  length  of  the 
interaction  region.  The  energy  resolution  obtained  in  this  mode  of  operation 
is  better  than  0.1  eV,  which  is  sufficient  for  studying  threshold  behavior 
and  fine  structure  on  the  ionization  curve.  At  higher  energies  where  no  fine 
structure  prevails  the  high  resolution  of  the  RPD  method  is  not  needed. 
Therefore  for  energies  above  30  eV  we  have  maintained  VR  constant  and 
used  the  ratio  I  /I  instead  of  Al./Al  in  Eq.  (38). 

"T*  — ’  T  “* 

The  potential  of  the  electron  collector  was  made  sufficiently  high 
to  ensure  saturation  of  electron  current  at  all  energies.  The  high  electric 
field  at  the  collector  surface  will  also  prevent  secondary  electrons  from 
entering  the  collision  region.  The  performance  of  the  electron  beam  is 
demonstrated  in  Fig.  14.  .  After  alignment  in  the  magnetic  field, 

the  current  is  independent  of  energy  to  within  +  1/2  %.  Typical  electron 
currents  used  in  this  investigation  were  from  0,1  -  1  pA  . 

The  ions  produced  between  the  parallel  plates  initially  have  zero 
kinetic  energy  and  will  move  to  the  upper  plates  in  a  proper  ExB  field 
configuration.  The  guard  plates  on  either  side  of  the  ion  collector  estab¬ 
lish  a  uniform  electric  field  E  and  accurately  define  the  length  &  of 
the  collision  region.  The  ion  trajectories  are  cycloidal,  with  drift  in 
the  direction  of  ExB  and  height 


E 


where  e  is  the  electronic  charge  and  M  and  Z  are  the  ionic  mass  and 
charge  number,  respectively.  As  long  as  half  of  the  parallel-plate  separation 
is  less  than  h  (i.e.,  ~  <  h  )  for  Z  =  Z  =3.  complete  collection  of  all 
ionic  species  is  assured.  This  requirement  was  always  satisfied  in  the 
experiment. 

The  actual  measurement  of  the  ion  current  I  is  complicated  by 
cesium  deposited  on  all  insulator  surfaces,  thereby  forming  a  leakage  path 
for  the  current.  By  using  ceramic  insulators  and  a  guard  ring  around  the 
ion  collector  feedthrough,  we  were  able  to  maintain  the  leakage  resistance 
of  this  part  above  10  Mfi  .  This  is  a  factor  of  100  -  1000  higher  than  the 
effective  input  impedance  of  the  Keithley  610  B  which  was  used  for  the  measure¬ 
ment  of  the  ion  current. 


4.2.3  Surface  Ionization  Detector  -  When  the  mean  free  path  of  the  atoms 
is  much  larger  than  characteristic  dimensions  of  the  apparatus,  Knudsen 
flow  can  be  assumed  and  the  density  n(T)  in  the  interaction  region  at  the 
temperature  T  can  be  found  from 


n(T)  =  N  - 


°'/W 


e'V 


(40) 


Here  N  is  the  gas  number  density  at  1  torr  at  Tq  =  273°K  and  plT^) 
is  the  vapor  pressure  at  the  surface  of  the  reservoir  which  has  temperature 

T_  .  Since  n(T)  «  T-'2  ,  it  is  more  than  sufficient  to  keep  the  temperature 

^  0 
T  in  the  collision  region  constant  to  within  +  5  K.  On  the  other  hand, 

the  strong  dependence  of  the  vapor  pressure  on  temperature  requires  that 

TCs  be  constant  within  +  0.05°K.  (A  temperature  change  AT^  =  0.1° 

corresponds  to  An/n  «  1%.) 


Experimentally,  prudence  must  be  exercised  in  the  application  of 

Eq.  (40).  First,  Knudsen,  or  molecular  flow  can  only  be  established  when 

27  28,29 

all  surfaces  are  covered  with  cesium.  Some  workers  have  reported 

that  this  might  take  many  days,  especially  at  surfaces  where  absorption  of 
atoms  will  obstruct  the  specular  reflection  of  atoms  that  is  required  in 
Knudsen's  theory.  Secondly,  there  exists  in  the  literature  disagreement  on 
the  values  of  tne  cesium  vapor  pressure  in  the  range  of  our  experiment. 
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By  comparison  we  find  that  the  recommended  values  of  Honig  and  of  Hultgren 

o  i  op 

et  al.  are  20  -  30%  below  Nesmeyanov's  average  values  and  Taylor  and 

25  33 

Langmuir's  formulas.  Recently,  Rozwadowsky  and  Lipworth  made  optical 

absorption  measurements  of  the  density  of  saturated  cesium  vapor  and  obtained 

34 

results  in  excellent  agreement  with  those  of  Taylor  and  Langmuir.  Con¬ 
sequently,  we  have  decided  to  use  the  latter  values  and  complement  these  with 
our  own  surface  ionization  detector  measurements. 

The  construction  of  the  surface  ionization  detector  (SID)  is  shown 
in  Fig.  15.  Above  the  center  of  a  guarded  circular  copper  plate  a  hot 
tungsten  wire  is  stretched  parallel  to  the  surface.  Positive  ions  produced 
over  the  length  D  are  collected  on  the  central  part  of  the  ion  collector 
and  give  a  current 

ISID=|n<o>F,  HI! 

where  (c)  is  the  average  thermal  velocity  of  the  cesium  atoms  and 
F  =  ndwD  .  By  comparison  with  another  SID  with  a  rectangular  collector, 
we  found  that  all  ions  generated  over  the  circular  collector  were  measured. 

The  tungsten  wire  was  aged  according  to  the  recommended  procedure 
24  o 

of  Taylor  and  Langmuir.  It  was  first  heated  to  2400  K  for  24  hours,  the 

temperature  was  then  increased  to  2600°K  for  a  period  of  1  hour,  and  finally 

the  wire  was  given  4-6  short  flashes  at  2900°K.  The  filament  temperatures 

were  determined  from  the  tables  of  Jones  and  Langmuir.  When  tungsten 

-9  -8 

filaments  were  treated  in  this  way  at  background  pressures  of  10  to  10 
torr,  we  found,  from  X-ray  diffraction  analysis,  a  layer  of  tungsten  carbide 
on  the  surface.  The  carbon  probably  comes  from  impurities  in  the  volume 
of  the  material.  Suppliers  of  tungsten  give  the  information  that  all  wires 
cotstain  some  carbon,  at  least  10  -  15  ppm,  because  of  the  fabrication  pro¬ 
cedure.  From  measurement  of  the  electron  emission  as  a  function  of  fila¬ 
ment  temperature  we  found  the  work  function  to  be  about  5.6  volts,  which  is 

sufficiently  high  to  ensure  ionization  of  cesium  atoms.  When  the  hi gh- 

-4  -5 

temperature  aging  took  place  at  much  higher  pressures  (10  -  10  torr),  the 

X-ray  diffraction  did  not  reveal  any  carbides  at  all.  The  explanation  of 
this  is  probably  that  carbon  at  the  surface  reacts  with  oxygen  to  give  CO. 
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15  Surface  ionization  detector.  The  diameter  of  the 


wolfram  wire  was 


d  =  1.28 


X  10”2  cm.  D  =  2.44  cm 


Nevertheless,  both  surfaces  with  and  without  carbides  gave  the  same  cesium 
densities  when  used  in  surface  ionization  detectors. 

Absolute  measurements  of  the  density  of  cesium  atoms  iri  the  vapor 

phase  by  means  of  surface  ionization  are  based  on  the  observation  of  Langmuir 

and  Kingdon^  that  all  atoms  leave  the  surface  as  ions.  (The  use  of  hot 

wires  for  density  determinations  in  atomic  beams  has  been  studied  and  dis- 
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cussed  by  Datz  and  Taylor  and  McFarland  and  Kinney.  )  By  using  Eq.  (41) 
for  determination  of  the  density  we  found  values  that  agree  with  those  cal- 
culated  from  Taylor  and  Langmuir's  formulas  to  within  +  Z%.  A  basic  require¬ 
ment  for  this  agreement  to  occur  is  that  the  whole  apparatus  be  in  thermo¬ 
dynamic  equilibrium  with  the  cesium  reservoir.  This  condition  was  reached 
about  two  weeks  after  cesium  was  introduced  into  the  system.  The  distance 
between  the  cesium  pool  and  the  top  of  the  tube  (about  1  meter) combined  with 
the  time  it  takes  to  cover  all  surfaces  with  cesium  accounts  for  part  of  the 
delay.  In  practice  we  found  that  the  density  increased  very  slowly  as  the 
equilibrium  condition  was  approached. 

4.2.4  Consistency  Checks  -  Some  of  the  precautions  that  must  be  made 

in  electron  beam  ionization  experiments  have  been  discussed  by  Briglia  and 
38 

Rapp.  First,  from  Eq.  (38)  we  would  expect  to  find  the  saturated  ion 
current  proportional  to  the  saturated  electron  current.  In  cesium  this  was 

verified  for  0.1  £  I  £  10  p A.  (A  test  case  in  helium  gave  I  «  I  for 

_10  -5  + 

10  £  I  £  10  A  .)  The  collection  efficiencies  for  electrons  and  ions 

have  already  been  discussed  in  Section  4.2.2. 

Furthermore,  proportionality  between  I  and  cesium  density  was 

10  12+  -3 

found  over  the  density  range  5  x  10  £  n  £  10  cm  .At  the  highest 

12  3  1  1^  15  1 

density  of  10  cm-0  the  electron  mean  free  path  is  A  =  (na)”  «  ( 10  “  x  10”  )~ 

3 

10  cm,  which  is  much  longer  than  the  length  l  of  the  interaction  region 
in  Fig.  13.  Therefore,  even  at  the  highest  density  used,  the  probability 
of  one  electron  producing  more  than  one  ionization  event  is  extremely  small. 

Beam  electrons  with  perpendicular  velocity  components  will  travel 
in  helices  parallel  to  the  magnetic  field.  The  path  length  along  one 
particular  helix  is,  of  course,  larger  than  the  length  A  as  defined  in 
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Fig.  13.  The  correction  formula  presented  by  Massey  and  Burhop  is  inappli- 

40 

cable  to  most  practical  guns,  as  pointed  out  by  Asundi.  The  electric  lens 
41 

configuration  in  our  electron  gun  tends  to  minimize  angular  spread  in  the 
beam.  Since  the  increase  in  path  depends  strongly  on  the  magnetic  field, 
the  best  way. to  check  the  constancy  of  the  ionization  length  is  to  measure 
I+/I  as  a  function  of  B  .  keeping  all  voltages  constant.  We  found  I  /I 
to  be  constant  for  the  range  of  parameters  used  in  the  experiment. 

Finally  we  checked  that  the  current  of  stray  electrons  to  the 
entrance  and  exit  electrodes  in  the  ionization  region  and  to  the  lower 
parallel  plate  was  1000  times  smaller  than  the  beam  current.  Stray  electrons 
were  not  detected  on  the  ion  collector. 

The  energy  scale  was  calibrated  by  letting  helium  into  the  tube 

42 

together  with  the  cesium  and  measuring  the  resonance  in  the  elastic  scatter 
ing  cross  section  at  19.3  eV.  The  resonance  shows  up  as  an  increase  in  trans 
mitted  electron  current.  The  other  point  on  the  energy  scale  was  taken  at 
the  appearance  of  Cs^  ions  at  3.89  eV.  In  addition,  we  plotted  the  electron 

90 

beam  current  as  a  function  of  energy  and  extrapolated  the  current  to  zero 
to  find  the  contact  potential  difference  between  the  collision  region  and 
the  retarding  electrode.  The  measured  voltage  difference  between  the  retard¬ 
ing  electrode  and  the  collision  region  includes  the  contact  potential  dif¬ 
ference  in  the  tube  as  well  as  contributions  from  the  leads  between  the  tube 
and  the  voltmeter.  This  total  correction  was  0.8  V,  independent  of  method 
of  measurement.  When  a  system  is  immersed  in  a  vapor,  all  surfaces  will  be 
automatically  covered  with  atoms  from  the  vapor,  thereby  providing  areas 
of  constant  contact  potential.  For  this  reason  we  believe  that  the  true 
contact  potential  difference  in  the  tube  is  much  smaller  than  the  measured 
value  of  0.8  V. 

The  spread  in  energy  in  the  electron  beam  can  be  attributed  to 
thermal  spread  as  the  electrons  leave  the  cathode  surface  and  to  non-uniform 
distribution  of  contact  potential  on  electrode  surfaces.  If  the  cathode 
temperature  is  Tc  ,  the  average  energy  of  the  emitted  electrons  is  2kT&  . 

The  Philips  impregnated  cathode  in  our  electron  gun  was  operated  at 
temperatures  of  about  800°K,  corresponding  to  an  average  thermal  energy 
of  0.14  eV.  We  obtained  information  about  the  spread  in  electron  energy 
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by  (1)  measuring  the  tailing  in  ion  current  and  by  (2)  retarding  potential 
measurement  of  the  current  to  the  electron  collector.  The  tailing  in  ion 
current  at  threshold  indicated  an  energy  spread  of  0.05  eV.  From  the  latter 
method  we  found  a  spread  in  energy  (defined  as  the  full  half-width  of  the 
differentiated  current-voltage  curve)  of  less  than  0.1  eV.  This  spread  was 
independent  of  the  beam  energy. 

4.3  Results  and  Discussion 

4.3.1  Data  Analysis  -  The  ion  current  that  is  measured  on  the  collector 
of  our  total  ionization  apparatus  is 

it  =  I(Cs+)  +  HCs**)  +  I(Cs+++)  +  ...  ,  (43) 

where  the  current  of  Csnf  ions  is 

I(Csn+)  =  I  nX0(Csn+)  .  (44) 

We  determine  the  current  of  singly  ionized  atoms  by  multiplying  the  total  ion 
current  with  the  ratio 

r  »  kcs-)  n\ot . 

This  ratio  is  identical  with  the  correction  factor  [Bio;  ]" 1  used  by  Brink, ^ 

n  J 

where 

a  =  a(Csn+) /a(Cs+)  . 
n 

18 

From  a  careful  analysis  of  Tate  and  Smith's  paper,  we  obtain  the  current 
ratio  shown  in  Fig.  16.  According  to  Tate  and  Smith  the  ionization  curve 
for  Cs  might  be  15%  too  low  due  to  incomplete  collection  of  ions.  It  should 
be  noted  that  the  Tate  and  Smith  ordinates  as  presented  in  their  paper  can 
be  interpreted  in  different  ways.  They  mention  expressions  such  as  "likelihood 
of  ionization",  "probability  of  ionization",  and  "efficiency  of  ionization", 
without  really  defining  any  of  these.  They  claim  the  ratio  of  the  different  ion 
currents  to  the  electron  current  gives  "a  measure  of  the  relative  probability 
of  the  production  of  the  different  ions."  As  will  be  shown  later,  different 
interpretations  will  give  different  results  for  the  Cs  cross  section. 
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FIG.  16  I(Cs+)/I^ot.  as  a  function  of  electron  energy 
as  obtained  from  Tate  and  Smith  (ref.  18). 
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4.3.2  Compilation  of  Results  -  Absolute  data  for  the  Cs+  cross  section  as 
a  function  of  electron  energy  are  shown  in  Fig.  17.  The  curve  in  the  middle 
represents  the  present  measurements.  Themeasurements  for  all  curves  were  per¬ 
formed  in  systems  in  which  cesium  vapor  was  in  supposed  thermodynamic  equil- 

22 

ibrium  with  a  liquid-phase  reservoir.  Korchevoi  and  Przonski  used  the 

43 

cylindrical  collector  arrangement  developed  by  Schulz,  whereas  Heil  and 
21 

Scott  modified  the  Tate  and  Smith-type  apparatus  by  depositing  the  ion 

collector  and  guard  ring  on  the  inside  of  a  glass  tube.  The  points  represent 

19  20 

crossed-beam  experiments  by  Brink  and  by  McFarland  and  Kinney.  If  McFarland 

and  Kinney  take  the  Tate  and  Smith  ordinates  as  probabilities  of  ionization. 

they  get  the  filled  circles  as  values  for  the  cross  section.  On  the  other  hand, 

if  the  Tate  and  Smith  ordinates  are  interpreted  as  ion  currents  proportional 

to  noICs”  ),  the  higher  values  indicated  by  the  open  circles  are  obtained. 

No  energy  values  have  been  assigned  to  the  maximum  cross  sections  reported  by 
19  20 

Brink  and  by  McFarland  and  Kinney.  From  the  shape  of  the  cross  section 
curves  we  see  that  a  rather  flat  maximum  occurs  at  28  eV.  Accordingly  the 
points  from  the  crossed-beam  experiments  have  been  plotted  at  this  energy. 

At  the  maximum  cross  section  (filled  circle  at  28  eV)  the  probability  of  ion- 
ization  is  approximately  equal  to  HCs  ) /I  ,  since  very  few  Cs  ions  are 
produced  at  this  energy.  The  McFarland  and  Kinney  value  at  maximum  is  there¬ 
fore  insensitive  to  interpretation  of  the  Tate  and  Smith  paper. 

Our  results  seem  to  coincide  with  those  of  Korchevoi  and  Przonski 
for  energies  up  to  about  10  eV.  Differences  in  initial  slope  are  not  dis¬ 
cernible  because  of  the  scale  used  in  Fig.  17.  (The  threshold  region  of  the 
cross  section  will  be  discussed  in  a  separate  paragraph.) 

It  is  interesting  to  normalize  the  curve  of  Tate  and  Smith  to  ours 
at  28  eV  and  compare  the  shape  of  the  cross  section  curves.  The  total  ion 

■j* 

current  I  as  reported  by  Tate  and  Smith  and  remeasured  here  falls  off 
almost  linearly  by  1/2  %  per  eV  from  28  to  100  eV.  Since  we  are  using  the 

-f  -J- 

ratio  of  HCs  from  Tate  and  Smith,  it  is  not  surprising  to  find 

good  agreement  in  this  energy  range.  Most  apparent  is  the  minimum  at  about 
19  eV  which  is  13%  below  the  maximum  on  the  Tate  and  Smith  curve  and  17% 
below  on  ours.  In  the  Heil  and  Scott  work,  the  "dip"  on  the  curve  is  about 
10%  of  the  maximum.  We  have  found  that  this  dip  decreases  with  increasing 
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ENERGY  (eV) 

FIG.  17  Ionization  cross  section  for  Cs+  vs  electron  energy.  The  fulldrawn 
curve  extending  to  100  eV  represents  the  present  measurements. 
Legend:  H+S:  Heil  and  Scott  (ref.  21). 

K  +  P:  Korchevoi  and  Przonski  (ref.  22). 

B:  Brink  (ref.  20). 

M  +  K:  McFarland  and  Kinney  (ref.  20). 

M:  McFarland  (ref.  45). 

T  +  S:  Relative  measurements  of  Tate  and  Smith  (ref.  18) 
normalized  to  the  present  maximum  at  28  eV. 
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electron  current,  presumably  because  higher  currents  result  in  a  larger  spread 
in  energy. 

Part  of  the  disagreement  between  the  various  results  can  be  under¬ 
stood  by  analyzing  the  experimental  techniques.  Since  Korchevoi  and  Przonski 
had  no  guard  rings  outside  their  ion  collector,  the  effective  length  of  the 
collision  region  is  larger  than  i>  .  The  cross  section  is  inversely  propor¬ 
tional  to  X  ,  which  might  explain  why  the  results  of  Korchevoi  and  Przonski 
are  higher  than  those  of  other  workers. 

The  lowest  cross  section  values  have  been  obtained  by  Heil  and  Scott. 

There  might  have  been  a  systematic  error  in  the  determination  of  the  cesium 

45 

density,  as  discussed  by  McFarland,  or  the  apparatus  might  not  have  been  in 
complete  thermodynamic  equilibrium. 

Basic  difficulties  in  crossed-beam  experiments  are  the  measurements 

of  the  cesium  atom  density  and  the  determination  of  the  interaction  volume. 

19  20 

Brink  and  McFarland  and  Kinney  were  using  two  different  experimental 
systems  with  different  geometry  factors.  In  spite  of  complicated  calibration 
procedures  the  results  agree  within  the  experimental  uncertainties,  which 
lends  confidence  to  the  reliability  of  this  type  of  experiment.  One  great 
advantage  of  the  crossed-beam  experiment  is  that  insulators  can  be  protected 
from  the  cesium  flux,  thereby  eliminating  the  leakage  current  problem  prevail¬ 
ing  in  vapor-phase  systems. 

The  accuracy  of  the  present  data  is  7 %,  based  upon  the  following 

estimates: 

AI+ 

-y1  =  5%;  £l_/I_  =  3%; 

An/n  =  3%;  and  A4/4  =  2%  . 

However,  the  high  energy  part  of  the  curve  might  be  too  low  due  to  a  system- 

-f.  *f 

atic  error  in  the  correction  factor  I(Cs  ) /I  .  The  maximum  systematic 
error  in  the  work  of  Tate  and  Smith  is  15%  but  this  error  might  be  energy 
dependent.  Since  McFarland  ^  and  Heil  and  Scott?*  have  been  using  the  same 
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correction  factor,  this  systematic  error  is  also  inherent  in  their  results. 

22 

Korchevoi  and  Przonski  avoided  this  ambiguity  by  limiting  their  measurements 
to  the  range  below  the  second  ionization  potential. 

By  taking  the  average  of  the  results  of  all  authors,  at  maximum 
°2 

we  obtain  a  value  of  9.8  A  .  It  has  been  assumed  here  that  the  curve  of 

°2 

Korchevoi  and  Przonski  extrapolates  to  12  A  at  28  eV. 

4-3.3  Ionization  Mechanisms 

The  shape  of  the  ionization  cross  section  curve  in  cesium  is  deter¬ 
mined  mainly  by  the  cross  section  for  removal  of  the  6s  electron.  According 
to  Geltman,  the  cross  section  close  to  threshold  (3.89  eV)  is  expected  to 
increase  in  proportion  to  the  electron  energy  in  excess  of  the  ionization 
energy.  In  electrical  discharges,  atoms  are  ionized  by  the  electrons  on  the 
"tail"  of  the  electron  energy  distribution,  and  for  this  reason  it  is  important 
to  know  the  initial  slope  of  the  cross  section.  Data  for  this  quantity  as 
well  as  the  range  of  linearity  are  shown  in  Table  V. 


Table  V.  Threshold  behavior  of  ionization  cross 

section. 

Author(s) 

00 

Slone  (A^/eV) 

Range 

above 

of  linearity 
threshold  (eV) 

Present  work 

2.7 

«  1 

Heil  and  Scott3 

2.2 

0.8 

Korchevoi  and  Przonski 

1.7 

unknown 

£ 

Nottingham 

5.7 

{v  1.5 

3  Ref.  21.  b  Ref.  22. 

C  Ref.  47. 

47 

The  measurement  of  Nottingham  was  made  in  an  electrical  discharge  where 

accurate  knowledge  of  the  parameters  needed  to  determine  the  cross  section 

is  difficult  to  establish.  It  is  surprising  to  observe  that  Korchevoi  and 

Pronzski,  who  reported  the  highest  maximum  cross  section  among  the  workers 

°2  , 

quoted  here,  obtained  the  lowest  value  of  the  slope  at  1.7  A  ,/eV.  Since  they 
do  not  state  the  range  of  linearity,  there  is  a  slight,  possibility  that  their 
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number  refers  to  the  tailing  part  of  the  ionization  curve  due  to  insufficient 

energy  resolution.  It  has  been  mentioned  already  that  the  maximum  of  our 

cross  section  agrees  closely  with  the  results  of  the  crossed-beam  experiments, 

°2 

and  for  this  reason  we  recommend  our  value  of  2.7  A  /eV  for  the  slope. 

As  a  first  attempt  to  explain  the  structure  (including  the  two 

maxima)  of  the  cross-section  curve,  we  have  calculated  the  partial  cross 

48 

sections  for  removal  of  6s,  5p,  and  5s  electrons  using  the  Gryzinsky  formula 

(Fig,  17).  The  total  cross  section  resulting  from  the  sum  of  the  partial 

Gryzinsky  cross  sections  agrees  very  closely  with  the  present  experimental 

results.  However,  the  most  conspicuous  difference  is  that  the  "dip"  on  the 

Gryzinsky  curve  is  only  1%  of  the  maximum,  whereas  a  value  of  17%  was  found 

here.  We  therefore  conclude  that  ionization  of  core  electrons  is  not  capable 

of  explaining  the  structure  of  the  ionization  curve,  and  that  additional 

mechanisms  are  needed.  Similar  types  of  structure  have  also  been  observed 

+  18  19  +  18  19  49  +  49 

on  the  ionization  curves  of  Hb  *  K  !  and  Mg  ,  but  in  the 

following  we  shall  restrict  our  discussion  to  Cs+  .  Heil  and  Scott  have 

suggested  that  the  first  peak  is  due  to  autoionization  and  the  second  due 

to  production  of  excited  ions.  The  Gryzinsky  calculation  shows  that  the 

°2 

cross  section  for  ionization  of  one  of  the  6p  electrons  is  about  0.5  A  , 
and  for  this  reason  the  production  of  excited  ions  seems  likely. 

02 

A  value  of  1  A  has  been  estimated  by  Heil  and  Scott  for  excitation 
of  the  autoionization  level.  Our  data  imply  a  value  of  about  1.6  K2  if  the 
first  peak  is  attributed  to  autoionization.  Interpreting  the  first  peak 
in  this  way  implies  the  existence  of  a  doublet  autoionization  level  (spin  = 
1/2)  at  about  12  eV  above  the  Cs‘  ground  state.  Feldman  and  Novick  have 
observed  a  quartet  (spin  =  3/2)  autoionization  level  at  12.6  +  0.3  eV. 

Their  state  is  metastable  (with  lifetime  about  40  nsec)  because  it  cannot 
make  a  dipole  transition  to  the  singlet  (spin  =  0)  continuum  states.  How¬ 
ever,  the  presence  of  a  quartet  state  at  this  energy  strongly  suggests  that 

5  2 

there  is  also  a  doublet  state  very  near  in  energy,  for  instance  the  6p  6s 

51 

electron  configuration  observed  by  Beutler  and  Guggenheimer.  Therefore, 
attributing  the  first  peak  to  autoinnization  is  probably  justified. 
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FIG.  18  Gryzinsky  calculation  of  the  cross  section.  Following  McFarland 

(ref.  46),  threshold  values  of  3.89,  17.2,  and  26  eV  were  used  for 
removal  of  6s,  5p,  and  5s  electrons,  respectively. 
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The  cross  section  for  excitation  of  the  doublet  autoionization 

level  is  expected  to  be  larger  than  for  excitation  of  the  quartet  metastable 

level  as  the  latter  requires  a  spin  change  of  the  5p  core  electron  of  the 

cesium  atom.  This  is  consistent  with  Feldman  and  Novick's  values  of  the  cross 
-2  -3  °2 

section  as  10  -  10  A  for  quartet  excitations  and  our  observation  of 

°2 

about  1.6  A  for  the  doublet  excitation.  The  lifetime  of  the  doublet  state 
is  too  short  to  be  observed  in  the  time-of-f light  type  experiment  of  Feldman 
and  Novick. 

In  conclusion,  it  should  be  pointed  out  that  the  ionization  mech¬ 
anisms  of  the  cesium  atom  are  not  completely  understood.  Experimentally,  one 
can  get  information  about  autoionization  by  studying  the  kinetic  energy 
(12.6  -  3.9  =  8.7  eV)  of  the  electron  ejected  from  the  autoionizing  level. 

This  will  be  the  objective  of  another  investigation. 
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Abstract 

Reported  electron  temperature  and  density  profiles  across  the  interelectrode 
space  of  a  thermionic  converter*  have  been  analyzed  u>ing  a  model  which  has  been  described 

o 

previously.  For  each  set  of  T  and  N  points,  effective  coefficients  for  ionization  and 

G  6 

for  radiative  energy  loss  are  computed.  These  coefficients  and  the  neutral  density 
(obtained  through  the  assumption  of  a  linear  heavy-particle  temperature  distribution 
across  the  gap,  and  the  ideal  gas  law)  yield  the  rates  for  net  ion  production  and  radia¬ 
tive  energy  loss  from  all  sources  except  the  resonance  lines.  Energy  loss  from  resonance 
lines  and  the  ion  production  cost  are  considered  separately.  The  following  general  con¬ 
clusions  can  be  drawn  from  the  consistent  behavior  of  the  results:  i)  the  area  within 
0.4  mm  of  the  emitter  is  a  source  of  ion-electron  pairs,  the  rest  of  the  plasma  being  a 
sink;  ii)  inelastic  energy  losses  decrease  monotonically  from  the  emitter  to  the  col¬ 
lector;  iii)  these  effects  are  relatively  independent  of  cesium  pressure,  but  are  generally 
proportional  to  output  current  density. 

Introduction 

Future  progress  in  the  effort  to  understand  the  operation  of  the  cesium  filled 
thermionic  converter  would  seem  to  depend  heavily  on  obtaining  an  accurate  solution  to 
the  equations  governing  the  transport  of  particles  and  energy  across  the  interelectrode 
plasma.  Two  of  the  most  important  terms  in  these  equations  are  the  rates  for  production 
or  loss  of  ion-electron  pairs,  and  the  rate  of  loss  of  electron  kinetic  energy  either 
directly  through  ionization  of  neutral  atoms  or  indirectly  through  escape  of  radiation 
from  the  volume  of  the  plasma. 

A  model  has  been  developed  which  yields  these  rates  if  the  following  informa¬ 
tion  is  given:  N  and  T  profiles  across  the  plasma,  electrode  temperatures  and  spacing, 
plasma  thickness,  and  cesium  pressure.  The  model  incorporates  the  assumptions  that  the 
plasma  is  radially  homogeneous;  that  the  only  species  present  are  ions,  electrons  and 
neutral  atoms;  that  mean-free-paths  for  electron-atom  collisions  are  small  compared  to 
plasma  dimensions;  that  gradients  in  T  and  N  are  negligible  in  distances  comparable  to 
these  mean-free-paths;  that  the  electron  velocity  distribution  is  Maxwellian;  and  finally 
that  the  plasma  is  optically  thin  to  all  but  resonance  radiation. 

The  model  consists  of  an  equation  for  the  rate  of  change  of  the  population  of 
each  atomic  level  due  to  inelastic  electron-ion  and  electron-atom  collisions,  and  spon¬ 
taneous  radiative  decay  and  recombination.  If  certain  criteria  are  satisfied,  the  excited 
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state  populations  can  be  assumed  constant  in  time,  and  the  rate  oi  change  ot  electron 

density  equals  the  rate  of  change  of  ground  state  density.  The  equations  can  then  be 

solved  for  the  excited  state  populations  in  terms  of  the  ground  state  density,  and  for 

the  rates  of  electron-ion  pair  production  and  radiative  energy  loss.  Since  the  model  has 

2 

already  been  discussed,  the  details  of  the  solution  will  not  be  considered  here.  An 

3 

exhaustive  presentation  will  become  available  shortly. 

The  Model 

Application  of  the  model  to  a  given  set  of  values  of  Ng  and  Tg  yields  results 
which  can  be  expressed  by  the  two  equations 

-dN  /dt  =  oN3  -  SN  N.  cnf3sec_1  (1) 

e  e  el 

-dR./dt  =  AN3  +  BN.  watt  cm"3  (2) 

1  e  1 

where  -dNg/dt  is  the  rate  of  loss  of  ion-electron  pairs, -dR^/dt  is  the  rate  of  loss  of 
radiative  energy  from  all  but  the  resonance  lines,  and  is  the  density  of  ground  state 
atoms.  The  coefficient  o'  can  be  interpreted  as  an  effective  three-body  collisional- 
radiative  recombination  coefficient,  and  S  can  be  interpreted  as  an  effective  ionization 
coefficient.  The  coefficients  A  and  B  do  not  have  such  a  simple  physical  interpreta¬ 
tion.  These  four  coefficients  are  functions  of  N  and  T  only. 

e  e  J 

It  should  be  clearly  understood  that  even  though  the  electron  density  at  any 
point  in  the  plasma  volume  is  constant,  the  rate  given  by  Eq.  (1)  is  not  necessarily  zero. 
This  rate  is  due  to  atomic  processes  only  and  must  be  balanced  by  macroscopic  processes 
such  as  diffusion.  It  should  also  be  noted  that  Eq.  (2)  represents  the  total  electron 
kinetic  energy  loss  only  when  Eq.  (1)  is  zero.  For  non-zero  -dNg/dt  there  is  an  addi¬ 
tional  energy  loss  term 

-  IQ 

-dR9/d'c  =  6.22  x  10  dN  /dt  watt  cm  (3) 

which  represents  a  cost  of  3.89  eV  per  ion-electron  pair  produced. 

Finally  (as  was  pointed  out  by  Kniazzeh  in  the  discussion  period),  while  it  may 
be  correct  to  assume  that  resonance  radiation  is  completely  trapped  in  solving  for  the 
coefficients  in  Eqs.  (1)  and  (2),  the  very  small  fraction  of  the  resonance  radiation 
which  does  escape  may  be  a  significant  energy  loss.  For  typical  conditions  the  Doppler 
and  pressure  contributions  to  the  line  half-width  are  comparable,  but  since  the  line 
center  is  completely  reabsorbed,  it  is  the  line  wings  which  determine  the  amount  of  radia¬ 
tion  escaping,  and  hence  only  the  pressure  broadening  need  be  considered/1  For  pressure 
broadened  lines  in  homogeneous  plasmas,  Holstein  gives  the  following  formulae  for  the 
fraction  of  radiation  (g)  which  escapes3 

gp  =  1.333(nkpL)"1/2  (4) 
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gc  =  1.376(TTkpR)  '  (5) 

for  plane  parallel  and  cylindrical  geometries  of  thickness  L  and  radius  R  ,  respectively. 
The  constant  is  the  absorption  coefficient  at  the  center  of  the  line  given  by^ 

1/2 

k  _  Jn_  /3a |  x-i  (6) 

kp  “  1.92  \  tuL  /  Xo 

The  statistical  weights  of  the  emitting  and  absorbing  state  are  and  u^,  respectively. 

For  the  first  resonance  line  of  cesium  (the  doublet  having  been  coalesced  into  a  single 

level),  Eqs.  (4)  and  (5)  give  g  =  2.074  x  10"3  L-^  and  g  =  2.141  x  10"3  R-1^, 

p  c 

where  L  and  R  are  in  cm.  Clearly,  for  spacings  greater  than  0.05  cm  more  than  99% 
of  the  resonance  radiation  is  trapped.  It  should  be  noted  that  this  result  is  dependent 
on  particle  pressures  or  temperatures  only  to  the  extent  that  pressure  broadening  must 
dominate  in  the  wings  of  the  line. 

7  -1 

The  transition  probability  of  the  6P  state  is  3.57  x  10  sec  ,  and  hence  an 
approximate  expression  for  the  power  loss  from  the  resonance  lines  in  a  cesium  plasma 
with  cylindrical  geometry  is 

-dRg/dt  =  1.72  x  10“14N2(6L_1/2  +  R"1/2)  watt  cm"3  (7) 

-3 

where  N2  is  the  6P  population  cm  and  €  is  the  emissivity  of  the  electrodes.  The 

contribution  from  higher  resonance  lines  can  safely  be  ignored  in  most  cases  because  of 

the  much  lower  density  of  these  levels. 

Solutions  of  the  set  of  rate  equations  yielding  the  rate  coefficients  in 

13  -2 

Eqs.  (1)  and  (2)  show  that  for  N  >  10  and  g9  <  10  ,  complete  resonance  trapping 

e  2 

can  be  assumed  without  introducing  significant  error.  Accordingly,  complete  sets  of  the 
coefficients  a,  S,  A,  and  B  were  computed  with  g  =  0.  These  results  were  seen  to  fit 
exponential  curves  of  the  type 

a  =  aQ  exp  [E^Cky 1/2]  (8a) 

S  =  SQ  exp  C_E2/kTe3  (8b) 

A  =  Aq  exp  [E3/(kTe)1/2]  (8c) 

B  =  Bq  exp  [-E4/kTg]  (8d) 


■j*  ^ 

The  numerical  constant  in  Eq.  (8)  of  Ref.  4  should  be  1.92  instead  of  2.74. 
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where  a  ,  etc.  and  E  ,  etc.  are  constants  which  are  dependent  only  on  N  .  Using 

0  i  V? 

Eqs.  (8)  and  interpolating  the  tabulated  values  of  a  ,  etc.  and  E^,  etc.  for  a  given 
N  reproduces  the  results  of  the  original  calculation  to  within  a  few  percent. 

The  degree  of  nonequi Ibrium  predicted  by  the  model  is  illustrated  in  Fig.  1, 
which  is  a  plot  of  the  results  of  Eq.  (1)  with  g  =  0  and  -dNg/dt  =  0.  The  ground  state 
density  predicted  by  the  model  is  plotted  vs  Ng  for  various  values  of  along  with 
that  predicted  by  the  Saha  equation 

N1  =  ( sr  r/2  Ne  ex» (l V'V  (9) 


where  E^  =  3.89  eV.  The  steady-state  nonequilibrium  calculation  is  seen  to  predict  a 
higher  value  of  for  a  given  N^,  reflecting  the  effect  of  radiative  decays  in  increas¬ 
ing  the  relative  population  of  the  lower  bound  levels.  A  more  conventional,  but  equiva¬ 
lent,  way  of  interpreting  these  plots  is  to  consider  that  the  model  predicts  a  lower  value 
of  N  for  a  given  and  Tg  due  to  depletion  of  high  bound  levels  and  the  continuum  by 
radiative  decay  and  recombination. 

Application  to  Experiment 

The  model  discussed  above  has  been  applied  to  experimentally  obtained  N  and 

,  1  6 
Tg  profiles  of  Reichelt.  These  data  were  taken  in  a  thermionic  diode  operating  in  the 

ignited  mode  with  an  interelectrode  spacing  of  1  mm  and  an  emitter-collector  radius  of 

8  mm.  The  emitter  and  collector  were  tungsten  and  nickel,  having  emissivities  of  ~  0.4 

and  Irrespectively,  at  their  operating  temperatures  to  radiation  in  the  neighborhood 

of  8500  A.  Cesium  pressures  were  above  1.9  torr.  Spectroscopic  measurements  of  Tg  and 

Ng  were  taken  at  four  points  across  the  gap,  each  separated  by  approximately  0.2  mm,  with 

a  spatial  resolution  of  ~0.1  mm.  Thus  each  point  is  distinct  from  the  adjacent  points 

and  over  a  volume  of  plasma  which  can  be  assumed  homogeneous.  The  measured  values  of 

Tg  and  Ng  ranged  from  2000  to  2800°K  and  from  5  x  10*3  to  5  x  lO^cm-3,  with  an  estimated 

experimental  error  of  5 %  and  10%,  respectively. 

At  these  operating  conditions,  the  assumptions  made  in  deriving  Eqs.  (8)  are 

—  _3  _3 

valid,  and  the  trapping  fractions  are  g  €  =  4.6  x  10  and  g  =  2.4  x  10  ,  implying 

p  c 

that  the  assumption  of  g  =  0  in  calculating  the  coefficients  in  Eqs.  (1)  and  (2)  is 

also  valid.  Due  to  the  high  degree  of  resonance  trapping,  which  forces  the  ground  and 

2 

first  excited  state  into  thermodynamic  equilibrium  at  T  ,  Eq.  (7)  can  be  simplified  to 
-dR3/dt  =  1.72  x  10"J'3  exp(-1.432/kT  )  watt  cm"3  (10) 


where  kTg  is  expressed  in  eV. 

The  experimental  values  of  Ng  and  Tg  were  used  to  generate  the  coefficients 
o',  S,  A  and  B  as  in  Eqs.  (8).  The  ground  state  density  was  computed  using  the  ideal 
gas  law 


P/k  =  N  T  +  N.T.  +  NT 
'  e  e  ii  a  a 


(11) 
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The  assumptions  were  made  that  the  ion  density  =  Ng  and  the  ion  temperature  T^  =  Tg, 
and  that  T  varies  linearly  between  the  electrode  temperatures.  The  total  atom  density 
N  is  related  to  by  a  simple  partition  function. 

The  results  obtained  for  -dNe/dt  are  shown  in  Figs.  3  and  4.  The  emitter  temper¬ 
atures  are  given  on  the  inset  J-V  curves.  Several  features  are  discernable:  1)  the  volume 
of  the  plasma  within  ~0.4  mm  of  the  emitter  appears  to  be  a  source  of  electron-ion  pairs, 
while  the  rest  of  the  plasma  is  a  sink;  2)  the  production  rate  near  the  emitter  is  gen¬ 
erally  independent  of  pressure  but  3)  is  generally  proportional  to  the  current  for  a 
given  emitter  temperature,  and,  4)  is  generally  inversely  proportional  to  emitter  temper¬ 
ature  for  a  given  current.  Departures  from  these  trends  reflect  departures  in  the  corres¬ 
ponding  relationships  between  Tg  and  output  current  J.  The  contribution  to  the 

electron  current  from  the  charged  particle  production  of  Figs.  3  and  4  is  negligible. 

20  -3-1 

A  production  rate  of  ~  5  x  10  cm  sec  across  the  entire  gap  would  be  required  to 
account  for  the  measured  current.  On  the  other  hand,  the  contribution  to  the  ion  current 
may  be  significant. 

These  conclusions  must  be  highly  qualified,  since  the  sensitivity  of  -dNg/dt  to 
T  is  great.  For  39  of  the  43  N  -T  sets  analyzed,  a  change  in  T  of  less  than  4%  is 
sufficient  to  yield  -dN  /dt  =  0.  For  the  remaining  points  the  required  change  in  T 

6  C 

was  less  than  10%.  Hence,  application  of  the  model  with  -dN  /dt  =  0  in  comparable  sys- 

v 

terns  can  be  expected  to  yield  electron  temperatures  within  experimental  error  for  known 

values  of  N  and  N  .  This  approach  has  been  used  by  Reichelt  and  others  with  Eq.  (1) 

X  6 

replaced  by  the  Saha  equation,  with  calculated  values  for  T  in  good  agreement  with 

e  13  -3 

experiment.  This  is  to  be  expected  from  the  plots  of  Fig.  1  for  N  >  5  x  10  cm  . 

C 

Use  of  the  Saha  equation  cannot,  however,  be  assumed  valid  at  lower  plasma  den¬ 
sities  where  departure  from  thermodynamic  equilibrium  is  greater.  The  sensitivity  of  the 

13  -3 

model  to  the  amount  of  resonance  trapping  increases  quickly  below  Ng  =10  cm  ,  so  that 
the  degree  of  nonequilibrium  can  be  expected  to  be  even  greater  than  that  shown  in  Fig.  1 
for  these  densities. 

The  results  obtained  for  -dR^/dt  are  shown  in  Figs.  4  and  5,  which  were  calcu¬ 
lated  using  the  experimental  values  for  Tg  .  The  energy  loss  due  to  radiation  from  all 
but  the  resonance  lines  is  seen  to  be  generally  proportional  to  current,  but  insensitive 
to  changes  in  emitter  temperature  or  pressure.  The  sensitivity  to  changes  in  Tg  is  not 
nearly  as  great  as  for  -dNft/dt  as. is  indicated  by  the  error  bars  on  one  set  of  data, repre¬ 
senting  changes  in  T  of  1%  and  5%.  The  integrated  total  loss  from  this  source  is  seen 
to  be  of  the  order  of  1  watt. 

Results  for  the  energy  loss  associated  with  net  ion  production  are  not  presented 

because  of  their  strong  dependence  on  T  .  It  should  be  noted,  however,  that  this  may 

®  18  —3  —1 

be  a  significant  source  of  power  loss,  since  a  production  rate  of  10  cm  sec  would  be 

-3 

at  a  cost  of  0.62  watt  cm 

Power  lo«'  including  all  radiation  is  plotted  in  Figs.  6  and  7.  These  data 
were  obtained  using  calculated  values  of  Te  t  i.e>,  under  the  assumption  that  -dNe/dt  =  0, 
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and  hence  does  not  include  any  contribution  from  ion-electron  pair  production  or  loss. 
Comparing  these  results  with  corresponding  Figs.  4  and  5,  it  is  seen  that  the  power  loss 
in  the  first  resonance  line  may  be  from  one-half  to  two-thirds  of  the  total  radiative  loss. 
The  results  presented  in  Figs.  6  and  7  can  be  compared  in  a  semi-quantitative 

7 

fashion  with  those  of  Nighan,  who  used  a  completely  different  approach  in  analyzing 

Q 

other  experimental  data  of  Reichelt.  The  conditions  of  emitter  temperature,  cesium 
pressure  and  output  current  were  similar,  but  the  emitter-collector  spacing  of  the  exper¬ 
iment  in  Ref.  8  was  2  mm.  Nighan  calculated  the  rate  of  electron  kinetic  energy  loss  to 

the  heavy  particles  due  to  both  inelastic  and  elastic  collisions,  obtaining  maximum  values 

-3  -2 

of  7.5  and  35  watt  cm  at  4.0  and  8.5  amps  cm  ,  respectively,  independent  of  cesium 
pressure. 
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Fig.  1  Ground  state  density  vs.  electron  temperature  Te  for  fixed  Ng. 
Nonequilibrium  (fte  =  0)  curves  ior  same  Ng  values  as  Saha  curves. 
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Fig.  2  Rate  of  electron-ion  pair  produc¬ 
tion  calculated  from  the  measured  Ne,Te 

profiles  and  the  nonequilibrium  model  vs  Fig.  3  Same  as  Fig.  2,  but  at  a 

distance  from  the  emitter.  The  collector  lower  cesium  pressure, 

is  at  1.0  mm.  The  J-V  curves  are  labeled 
with  the  emitter  temperature. 
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Fig.  4  Power  loss  from  all  sources  except  the  resonance 
from  the  measured  Ne,Te  profiles  and  the  nonequilibrium 
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Fig.  5  Same  as  Fig.  4,  but  at  a  lower  cesium  pressure. 

indicate  the  effect  of  charges  in  T  of  1%  and  5% 
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Fig.  6  Total  radiative  power  loss  calculated  from  the  measured  N 
profiles  and  the  nonequilibrium  model  with  N  =0,  R  =R. +RQ  .  e 
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Fig.  7  Same  as  Fig.  6,  but  at  a  lower  cesium  pressure. 
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